Optimal Radio Window for the Detection of Ultra-High-Energy Cosmic Rays and 

Neutrinos off the Moon. 
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When high-energy cosmic rays impinge on a dense dielectric medium, radio waves are produced 
through the Askaryan effect. We show that at wavelengths comparable to the length of the shower 
produced by an Ultra-High Energy cosmic ray or neutrino, radio signals are an extremely efficient 
way to detect these particles. Through an example it is shown that this new approach offers, for 
the first time, the realistic possibility of measuring UHE neutrino fiuxes below the Waxman-Bahcall 
limit. It is shown that in only one month of observing with the upcoming LOFAR radio telescope, 
cosmic-ray events can be measured beyond the GZK-limit, at a sensitivity level of two orders of 
magnitude below the extrapolated values. 
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I. INTRODUCTION 

The interest in determining the flux of Ultra-High En- 
ergy (UHE) cosmic rays and neutrinos is manyfold. The 
origin of the highest energy cosmic rays is a major re- 
search topic as the existence of these particles requires 
very spectacular events on a cosmic scale. At energies 
beyond the so-called Greisen-Zatsepin-Kuzmin (GZK)- 
limit Q, (at an energy of about 6 lO^^ eV) the spec- 
trum of cosmic rays is expected to drop rather drasti- 
cally. The mechanism for this is that at these energies 
the cosmic rays produce pions when scattering off the 
microvifave background while traversing distances of the 
order of 10 Mpc. The existence of this cutoff in the spec- 
trum has not been verified unambiguously up to now. 
The flux of cosmic rays beyond the GZK cutoff deter- 
mines the sources for UHE cosmic rays within a range 
of about 10 Mpc H i.e. close on astronomical scales. 
This is a very exciting prospect as presently there are no 
known sources at this proximity. The presently proposed 
method is very efficient for determining this fiux. 

Another point of interest lies in the detection of UHE 
neutrinos. These could be created by UHE protons pro- 
ducing 7r+ mesons when scattering off the microwave 
background (the GZK mechanism as mentioned above) 
which, through weak decay, produce neutrinos. These 
GZK neutrinos have thus far never been observed. 

There are also other, more speculative, models predict- 
ing UHE neutrinos. These models belong to a generic 
class known as top-down (TD) models, where UHE par- 
ticles owe their origin to the decay of some supermassive 
X-particle of mass mx- Their decay products, the UHE- 
cosmic rays, can have energies up to mx- These mas- 
sive X particles could be topological defects or magnetic 
monopoles that could be produced in the early Universe 
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during symmetry-breaking phase transitions envisaged in 
grand unified theories (GUTs); see [ESS Hi for reviews. 

As an efficient method to determine the fluxes of UHE 
particles we are investigating the production of radio 
waves in a particular frequency window when a UHE 
particle hits the moon. Askaryan predicted as early as 
1962 Q that particle showers in dense media produce co- 
herent pulses of microwave Cerenkov radiation. Recently 
this prediction was confirmed in experiments at accel- 
erators |lOl | and extensive calculations have been per- 
formed on the development of showers in dense media 
to yield quantitative predictions for this effect The 
Askaryan mechanism lies at the basis of several exper- 
iments to detect (UHE) neutrinos using the Cerenkov 
radiation emitted in ice caps [H m, salt layers [l|, 
and the lunar rcgolith. The pulses from the latter pro- 
cess are detectable at Earth with radio telescopes, an idea 
first proposed by Dagkesamanskii and Zheleznyk |l5l | and 
later by others . Several experiments have since been 
performed llSj to find evidence for UHE neutrinos. 
All of these experiments have looked for this coherent 
radiation near the frequency where the intensity of the 
emitted radio waves is expected to reach its maximum. 
Since the typical lateral size of a shower is of the order 
of 10 cm the peak frequency is of the order of 3 GHz. 

Here we propose a different strategy to look for the 
radio waves at considerably lower frequencies where the 
wavelength of the radiation is comparable in magnitude 
to typical longitudinal size of showers. It has been noted 
before [18] that a new generation of low-frequency digital 
radio telescopes will provide excellent detection capabil- 
ities for high-energy particles, thus making our consider- 
ation here very timely. We show that the lower intensity 
of the emitted radiation, which implies a loss in detection 
efficiency, is compensated by the increase in detection ef- 
ficiency due to the near isotropic emission of coherent 
radiation. The net effect is an increased sensitivity by 
several orders of magnitude, for the detection of UHE 
cosmic rays and neutrinos at frequencies which are one 



2 



or two orders of magnitude below that where the inten- 
sity reaches its maximum. At lower frequencies the lu- 
nar rcgolith becomes increasingly transparent for radio 
waves. This implies for the detection of UHE neutri- 
nos that there are two gain factors when going to lower 
energies; i) Increased transparency of the lunar regolith 
already stressed in Ref. jl^ , and ii) Increased angular ac- 
ceptance, stressed in this work, which gives much larger 
count rates. 

In the following two sections we explain quantitatively 
the idea of an optimum frequency by taking cosmic-ray 
and neutrino-induced radio-emission from the Moon as 
a specific example. The advantage of going to lower fre- 
quencies also applies to other experiments where the radi- 
ation crosses a boundary between a dense medium to one 
with a considerably lower index of refraction. In Section 
IV we propose two specific observations, one for an ex- 
isting facility, the Westerbork Synthesis Radio- Telescope 
array (WSRT), and one for a facility which will be avail- 
able in the near future, the Low-Frequency Array (LO- 
FAR). 



II. MODEL FOR RADIO EMISSION 

There exist two rather different mechanisms for radio 
emission from showers triggered by UHE cosmic rays or 
neutrinos, where each has received considerable atten- 
tion recently. One is the emission of radio waves from 
a shower in the terrestial atmosphere. Here the pri- 
mary mechanism is the synchrotron acceleration of the 
electrons and positrons in the shower due to the geo- 
magnetic field, called geosynchrotron radiation, which 
has recently been confirmed with new digital radio tech- 
niques 0, 123, 1^ [12, ■ The second mechanism ap- 
plies to showers in dense media, such as ice, salt, and 
lunar rcgolith, where the front end of the shower has a 
surplus of electrons. Since this cloud of negative charge 
is moving with a velocity which exceeds the velocity of 
light in the medium, Cerenkov radiation is emitted. For a 
wavelength of the same order of magnitude as the typical 
size of this cloud, which is in the radio-frequency range, 
coherence builds up and the intensity of the emitted ra- 
diation reaches a maximum. This process, known as the 
Askaryan effect is the subject of this paper. 

The intensity of radio emission (expressed in units 
of Jansky's where 1 Jy = 10^^^ W m^^Hz^^) from a 
hadronic shower, with energy Es, in the lunar regolith, 
in a bandwidth Av at a frequency v and an angle 0, can 
be parameterized as (see Appendix 
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where I/O = 2.5 GHz [3, d is the distance to the observer, 
and dmoon = 3.844 x 10® m is the average Earth-Moon 
distance. The angle at which the intensity of the radia- 
tion reaches a maximum, the Cerenkov angle, is related to 
the index of refraction (n) of the medium, cos 0c = ^/n. 
Crucial for our present discussion is the spreading of the 
radiated intensity around the Cerenkov angle, given by 
Ac (in degrees). The sin0 factor in Eq. reflects the 
projection of the velocity of the charges in the shower 
on the polarization direction of the emitted Cerenkov ra- 
diation. The dependence of Z as defined in Eq. Q is 
suggested by working out some specific cases i see Ap- 
pendix^ For small values of Ac it coincides with the for- 
mula Eq. ljAl|l found in much of the hterature [THIT^IT^ 
however, Eq. |^ is more accurate for large spreading an- 
gles. 

The spreading of the radiated intensity around the 
Cerenkov angle, Ac, is, on the basis of general physical 
arguments, inversely proportional to the shower length 
and the frequency of the emitted radiation. Based on 
the results given in Ref. it can be parameterized as 
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where L{Es) is the shower length which depends on the 
energy. In Ref. [l^l, calculated results are given for the 
shower length (in units of radiation lengths, equal to 
22.1 g/cm^ ^3 for lunar regolith) which at the highest 
energies can be parameterized as 



L{x) = 12.7 



(4) 



where x = logjQ(£'/102° eV). At an energy of 10^° eV 
this corresponds to a shower length of approximately 
1.7 m, where we take the density of the regolith to be 
approximately 1.7 g/cvc?. For a frequency v = 200 MHz, 
where the wavelength is of the same order as the shower 
length, we should expect on general arguments that the 
radiation spreads over an angular range (2Ac) that is 
comparable to the Cerenkov angle, 9c = 56°. This in- 
deed corresponds to the value Ac = 21.52°, obtained 
from Eq. In Appendix 1X1 the predictions for the an- 
gular spread based on the parametrization Eq. ^ using 
Eq. is compared with analytical calculations for some 
(simplified) shower profiles showing excellent consistency 
with the calculated shower length. 

In our simulations we have taken into account the at- 
tenuation of radio waves in the regolith. As a mean value 
for the attenuation length for the radiated power we have 
taken = (9/t/[GHz]) m |2^, which is the same value 
as used in the analysis of the GLUE experiment [l8l |. 
This value is obtained from loss-tangent measurements 
performed on samples of lunar basalt brought back from 
the Moon . The measured values show a rather large 
variation, the effects of which are investigated in Sec- 
tion In principle the layer of regolith is only 10-20 m 
thick under which there a thick layer of fractured rock 
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going over into solid bedrock. As shown in Appendix 
the bedrock is as efhcient in emitting radio waves at the 
lower frequencies as the regolith. At higher frequencies, 
due to the larger attenuation length, only the relatively 
thin upper layer contributes to radio emission. All in 
all this implies that for the calculation of the acceptance 
the structure of the deeper layers (rock v.s. regolith) is 
not important. In the calculations we have therefore in- 
cluded radiation coming from a depth of at most 500 m 
treating for simplicity, and without loss of accuracy, the 
whole layer as behaving like regolith. It is argued in Ap- 
pendix [ni that this indeed gives a realistic estimate for 
the acceptance calculations. 

In the calculations for cosmic-ray-induced showers we 
assumed that the shower occurs effectively at the lunar 
surface. As argued in Appendix ^ only a very small 
depth is necessary for Cerenkov radiation to be emitted. 
For neutrino-induced showers an energy-dependent mean 
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free path jH has been used, = 130 ( ^"g,^*^^ ) km 
which is appropriate for regolith. 

A crucial point in the simulation is the refraction of 
radio waves at the lunar surface as was already stressed 
in Ref. Since the index of refraction of the lunar 

regolith is relatively large, n = 1.8 corresponding to a 
Cerenkov angle of 6*^ = 56°, much of the radio wave 
which is emitted by the shower suffers internal reflection 
at the surface. Only radiation with an angle of 90° — Oc 
or less with respect to the normal to the surface will be 
emitted from the Moon. Since most showers, being cos- 
mic ray or neutrino induced, are directed towards the 
center of the Moon, internal reflection will severely di- 
minish the emitted radiation at high frequencies where 
the Cerenkov cone is rather narrow. The major advan- 
tage of going to lower frequencies is that the spreading 
Ac around dc increases, allowing for the radiation to es- 
cape from the lunar surface. With decreasing frequency 
the peak intensity of the emitted radiation decreases (see 
Eq. (jAll) ). however the peak intensity increases with in- 
creased particle energy. The net effect is that at suffi- 
ciently high shower energies the aforementioned effect of 
increased spreading is far more important, resulting in a 
strong increase in the detection probability. 

To be able to address these issues quantitatively we in- 
troduce the detection efficiency 'D{E,v). It is defined as 
the probability that a cosmic ray (or neutrino) hitting the 
Moon (at an arbitrary angle and position) with energy 
E would produce radio waves at frequency v which is de- 
tectable at Earth. We regard an event to be detectable 
when the power of the signal is 25 times larger than the 
"noise level" which we take equal to Fnoise ~ 20 Jy [2^, 
i.e. a minimal detection level of 500 Jy is taken using a 
bandwidth of Ai/ = 20 MHz in Eq. |(T}. These values 
can be considered typical for LOFAR, see Section IIVI 
It is illustrative to plot the differential detection proba- 
bility, dV/dfl, for a particular area dA = R^dfl on the 
Moon where R is the lunar radius. In Fig. ^ dV/dil is 
plotted for cosmic rays at an energy of 4 x 10^^ cV for 
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FIG. 1; Differential detection probability, dV/dil, for a cos- 
mic ray of energy 4 x 10^^ eV hitting the Moon as function of 
apparent distance from the center of the Moon, b, for different 
detection frequencies. 



TABLE I: Calculated detection probabilities for cosmic-ray 
induced showers at Ecr — 4x 10^^ eV, labelled as 2?cr(40), and 
for neutrino-induced showers at Ev = 2 x 10^^ eV, labelled as 
I'^(200), for different frequencies u (in units of GHz). 
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different frequencies as a function of the relative distance 
from the center of the face of the Moon, b, where 6=1 
corresponds to the rim of the Moon. At lower frequen- 
cies the length of the shower becomes comparable to the 
wavelength of the radiation. Coherent emission of radio 
waves thus happens over a large angular range instead of 
only within a narrow cone at the Cerenkov angle. The 
effect of this is that even cosmic rays hitting the center 
of the face of the Moon (as seen by us) at very oblique 
angles are detectable. If the radiation were emitted in a 
very narrow cone around the Cerenkov angle, much of the 
radiation, even if the cone is directed towards the Earth, 
would be be internally reflected off the lunar surface. All 
this implies that at high frequencies only the rim of the 
Moon contributes to dV/dfl and that this contribution is 
not really large since only very inclined cosmic rays may 
produce detectable emission. At lower frequencies the 
whole surface of the Moon contributes with a relatively 
large probability since a large range of angles contribute. 
This trend is clearly discernable in Fig. ^ 

In Table m the detection probabilities, D {dV/dil in- 
tegrated over the lunar surface), are given as a func- 
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tion of frequency. This shows a strong increase of T> 
with decreasing frequency. For a given shower the radi- 
ation that is transmitted through the surface (which is 
not internally reflected) is proportional to , where the 
quadratic dependence is due to the fact that the phase 
space is in both the polar and in the azimuthal angle. 
An additional factor of Ac comes from the fact that also 
for showers making an angle of up to Ac with respect 
to the tangent to the surface, radio waves may be trans- 
mitted through the surface. In total one expects thus 
T) oc Ac oc which agrees rather well with the num- 
bers given in Tabled 

For neutrino-induced showers only 20% of the initial 
energy is converted to a hadronic shower, while the re- 
maining 80% is carried off by the Icpton. This energetic 
lepton will not induce a detectable radio shower. For a 
muon, the density of charged particles will be too small, 
while the shower of a UHE electron will be extremely 
elongated due to the Landau-Pomcranchuk-Migdal ef- 
fect [23. The width of the Cercnkov cone will thus 
be very small which makes the shower practically unde- 
tectable. For the present calculations we therefore have 
limited ourselves to the hadronic part of the shower which 
carries 20% of the energy of the original neutrino. To be 
able to compare the results for cosmic-ray and neutrino- 
induced showers in Table|l]the latter have been calculated 
for a 5-timcs higher energy of the incoming particle. 

While the showers for hadronic cosmic rays develop 
very close to the lunar surface, those of UHE neutrinos 
are distributed over a rather large range of depths from 
the surface. As a result of the long attenuation length 
for neutrinos in matter (tens of km at our energies), a 
large fraction of the neutrinos create showers too deep 
inside the Moon so that the radio waves are attenuated 
to below the detectable threshold at Earth. Roughly, 
one thus expects that the T) for neutrino- induced showers 
will be a factor Xr/Xu smaller than V for cosmic-ray- 
induced showers at the same shower energy. This factor 
explains much of the difference between 'Dcr and 2?i/ given 
in Table n 

A deviation from the simple Xr/Xi, scaling occurs since, 
due to the relatively long neutrino mean free path, the 
hadronic part of the neutrino-induced shower has a cer- 
tain probability of being directed towards the surface 
(in contrast to showers induced by hadronic cosmic rays 
which arc always directed into the Moon) . This increases 
the probability that part of the radio signal is emitted 
from the Moon. The latter can be seen from Fig. |2] where 
at 2.2 GHz the acceptance ring is broader than in Fig.^ 

An additional advantage of using lower frequencies is 
that the sensitivity of the model simulations to large- 
or small-scale surface roughness is diminished. Since at 
lower frequencies already a sizable fraction of the radia- 
tion penetrates the surface, its roughness will not make a 
major difference. This is in contrast to high frequencies 
where most of the radiation is internally reflected when 
surface roughness is ignored. 

The increased spreading of the radiation around the 
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FIG. 2: Same as Fig.Qbut for showers induced by a neutrino 
of energy 2 x lO'^^ eV. It is assumed that 20% of the neutrino 
energy is deposited in an hadronic shower, the shower induced 
by the high-energy lepton is ignored in the present calculation 
since its Cerenkov cone is very sharp and thus a negligible 
fraction of its energy will penetrate through the lunar surface. 



Cerenkov cone at lower frequencies strongly increases the 
detection efficiency but at the same time decreases the 
sensitivity to the direction of the original cosmic ray or 
neutrino. Part of this can be recovered by measuring the 
polarization direction of the radio waves. The electric 
field is 100% polarized in the direction of the shower, 
which has been confirmed in laboratory experiments |lOl | . 



III. DETECTION LIMITS 

In Fig. 121 the detection limits for UHE cosmic rays for 
different radio-frequency ranges arc compared with data 
from the AGASA [2^ (points) and a linear extrapola- 
tion based on the data from the HiRes experiment 
(grey bar). The model-independent limit is defined as 
dN/dEiim = Q/{V X h) where h equals the observation 
time in hours and Q = 1.16 x 10~^^ cm~^s~^sr~^ is the 
full phase-space for 1 hour on the Moon, assuming an 
isotropic distribution for the cosmic rays or neutrinos and 
assuming the whole face of the Moon is in the antenna 
field of view. As in the previous section we have calcu- 
lated the detection probability T) for a signal threshold 
of 500 Jy at all frequencies. It should be noted that for 
v ~ MHz there is a strong increase in the sky temper- 
ature and we have used a ten-fold higher threshold. As a 
result of the higher detection threshold the flux limit lies 
considerably higher. 

From Fig. O one clearly sees that with decreasing fre- 
quency one loses sensitivity for lower-energy particles. 
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FIG. 3: Flux limits (assuming a null observation) for cos- 
mic rays as can be determined in a 100 hour observation (see 
text). In the curves for u = 30 MHz a ten fold higher detec- 
tion threshold is used, corresponding to the higher sky tem- 
perature at this frequency. The points given correspond to 
the data of the AG AS A experime nt [231 . the grey band is an 
extrapolation of the HiRes data l2g|[ . The thick black line 
corresponds to the best possible limit (vanishing detection 
threshold) . 



This follows directly from Eq. (pj since with decreasing 
frequency the maximum signal strength decreases and 
thus one exceeds the detection threshold only for more 
energetic particles. If the energy of the cosmic ray is 
more than a factor 4 above this threshold value the anal- 
ysis presented in the previous section applies and the 
detection limit improves rapidly with decreasing radio- 
frequency until one reaches a frequency of 100 MHz where 
one obtains the optimum sensitivity. Decreasing the 
frequency even lower provides no gain since the detec- 
tion limit has already reached the optimum dN/ dEopt = 
(5/(0.5 X ft.), given by the heavy black line in Fig.|3| 

In Fig. 0] we compare the detection limits for UHE neu- 
trinos at different frequencies with the results obtained 
from the GLUE experiment . One sees similar trends 
as in the predictions for cosmic rays, in particular the 
large gain in the determined flux limits with decreasing 
frequency. At higher energies the limits for neutrinos do 
not increase as steeply as those for cosmic rays. This is 
because the neutrino mean-free-path decreases with en- 
ergy, therefore increasing the probability for the neutrino 
to initiate a shower close to the surface where the atten- 
uation of the radio waves is small. Our result at 2.2 GHz 
happens to lie close to that of the GLUE experiment. 
The dashed-dotted curve in Fig. ^ shows the results of 
a calculation where we have reproduced the simulation 
for the GLUE experiment, i.e. included (in a somewhat 



FIG. 4: Similar flux limits as shown in Fig. Elbut for UHE 
neutrinos. The open squares are the limits determined from 
the GLUE experiment p^ . 

simplified manner) the effects of averaging over lunar- 
surface slopes of 10°, including a 10% coverage of the 
Moon, and have used the appropriate detection thresh- 
old. This result lies close to the published limits of the 
GLUE experiment. 

In this work we mainly address the detection limits for 
fluxes of UHE cosmic rays or neutrinos and a distinction 
between the two does not have to be made. When one 
would measure one or more events the question of dis- 
tinguishing between the two kinds becomes interesting. 
For a single event the frequency dependence of the pulse 
(assuming a broad band acceptance) might give an indi- 
cation since for deep showers the high-frequency part will 
suffer a larger attenuation than the low- frequency part. 
For a large number of events the distribution over the 
lunar surface could be exploited. 



IV. WSRT &j LOFAR PREDICTIONS 

The Westerbork Synthesis Radio Telescope (WSRT) 
consists of fourteen 25 m parabolic dishes located on an 
east- west baseline extending over 2.7 km It is nor- 
mally used for super-synthesis mapping, but elements of 
the array can also be coherently added to provide a re- 
sponse equivalent to that of a single 94 m dish. Observing 
can be done in frequency bands which range from about 
115 to 8600 MHz, with bandwidths of up to 160 MHz. 
The low frequency band which concerns us here covers 
115-170 MHz m . Each WSRT element has two receivers 
with orthogonal dipoles enabling measurement of all four 
Stokes parameters. In tied-array mode the system noise 
at low frequencies is Fnoise =600 Jy. To observe radio 
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FIG. 5: Flux limits on UHE cosmic rays as can be determined 
in a 30 day observation with the LOFAR antenna system and 
a 100 hour observation with WSRT. The data shown are the 
same as in Fig.|3 



bursts of short duration, the new pulsar baekend (PuMa 
II) will be used. It ean provide dual-polarization base- 
band sampling of eight 20 MHz bands, enabling a maxi- 
mum time resolution of approximately the inverse of the 
bandwidth. In the configuration which we propose to 
use, four frequency bands will observe the same part of 
the moon with the remaining four a different section. In 
total, coverage of about 50% of the lunar disk can be 
achieved. 

An even more powerful telescope to study radio flashes 
from the moon will be the LOFAR array ■ With a 
collecting area of about 0.05 km^ in the core (which can 
cover the full moon with an array of beams), LOFAR will 
have a sensitivity about 25 times better than that of the 
WSRT. LOFAR will operate in the frequency bands from 
30-80 and 115-240 MHz where it will have a sensitivity 
of about Fnoise =600 Jy and Fnoise =20 Jy, respectively. 
The Galactic background noise will become the dominant 
source of thermal noise fluctuations at frequencies below 
about 100 MHz. It therefore appears that the optimal 
radio window for the detection of cosmic ray or neutrino 
induced radio flashes from the moon will be around 100- 
150 MHz. 

The probability for a observing a signal simultaneously 
in four frequency bands with a power of 20x Fnoise per 
frequency band, where Fnoise is the noise level per po- 
larization direction, as a random fluctuation of the back- 
ground noise is less than 0.001 (equal to 3a signiflcance) 
for a 100 h observing period at a sampling rate of 40 
MHz, typical for our PUMA-II back-end. Simulations 
show that a pulse of intensity 25xFnoise: interfering with 
the noise background, can be detected with 3cr signifi- 



cance at a probability greater than 80%. For this reason 
we have assumed in the calculations a detection thresh- 
old of 25xFnorse for both the WRST and the LOFAR 
telescopes. 

A simulation for LOFAR, taking ly = 120 MHz, band- 
width of Ai^ ~ 20 MHz, a signal-detection threshold of 
500 Jy, and an observation time of 30 days is shown in 
Fig. [SI for cosmic rays and in Fig. for neutrinos. The 
results are compared with the limits that can be obtained 
from a presently proposed observation for 100 hours at 
the WSRT observatory assuming a detection threshold 
of 15,000 Jy, V = 140 MHz, bandwidth of Ai/ = 20 MHz, 
and a 50% Moon coverage. 

In Fig. 13 the predicted spectrum of GZK neutrinos is 
taken from Ref. [sj. The prediction of top-down (TD) 
neutrinos is based on the decay of a topological defect 
with a mass of 10^"' eV as was calculated in Ref. |^ l3^. 
The Waxman-Bahcafl (WB) hmit |10|, 

based on theo- 
retical arguments, is an upper limit on the neutrino flux 
which is consistent with the data on the fluxes of UHE 
cosmic rays. The limits are also compared with those 
from the GLUE [l3| and FORTE experiments which 
are calculated as model independent limits similar to our 
limits. The limit from the RICE [s^l experiment has been 
calculated assuming different power-law spectra for the 
neutrinos. In general such a limit lies below the model- 
independent limit [l^ . 

With the existing WSRT a limit on the neutrino flux 
can be set which falls just above the WB bound. How- 
ever, even this will constrain different top-down scenar- 
ios, discussed in the literature. With the proposed LO- 
FAR facility this limit can be improved considerably to 
reach, for the first time, a limit well below the WB bound 
for neutrinos. In addition one has a good chance to see 
evidence (in only a 30 day period) of cosmic ray events at 
an energy one order of magnitude higher than presently 
observed. 



V. ACCURACY OF THE PREDICTIONS 

In this section we address the robustness of our ac- 
ceptance calculations. Some of the sources of model- 
dependence in the calculation have already been dis- 
cussed. A potentially large one is due to the uncertainty 
in the thickness of the regolith layer. As is argued in 
Appendix [nj the final result is rather insensitive to this 
thickness since the rock below the regolith is as efficient, 
if not more so, in emitting radio- waves at the frequencies 
we arc interested in. 

The loss tangent of the regolith, which determines the 
attenuation of radio waves, is found to be dependent 
on the metallic composition of the regolith. Different 
Apollo samples show a large variation. The extremes for 
radio-attenuation distance in regolith vary between 2/1/ 
m/GHz to about 25/i^ m/GHz |2J|. We have used an in- 
termediate value of m/GHz for our estimates. As can 
be seen from Fig. [7| (please note the expanded scale) the 
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FIG. 6: Flux limits on UHE neutrinos as can be determined 
with WSRT and LOFAR observations (see text) are com- 
pared with various models, in particular, WB |30l (vertical 
bars), GZK IH (dotted thin line), and TP |l l32 | (solid thin 
line). Limits from the RICE GLUE [H, ANITA 
and FORTE experiments are also shown. 
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for WSRT. Note the different scale in this figure. 



two extremes for the radio-absorbtion distance result in 
roughly a variation in the acceptance for neutrinos which 
is equal to the variation in the radio-attenuation length. 
The reason is that the thickness of the layer of the lu- 
nar crust which can be 'seen' on earth is proportional to 
the attenuation length of radio waves. It should be real- 



FIG. 8; Model dependence of flux limits on UHE cosmic rays 
for WSRT (and the LOFAR antenna system). 



ized that the showers initiated by hadronic cosmic rays 
are close to the surface for which the exact value of the 
loss-tangent is unimportant. 

We show in Appendix^that the angular spread of the 
emitted radio- waves from a shower depends on the length 
of the shower. This length is proportional to the stopping 
power Xq. As the stopping power is mostly determined 
by electronic processes it is strongly dependent on the el- 
ements in the rock or regolith. We have determined the 
stoppin g po wer for some of the Apollo samples analyzed 
in Ref. |43 . li^ l . This shows that the actual variation 
in the stopping power is not as large as one might have 
expected. We have determined Xq = 22.1 g/cm^ for the 
'A-17 HIGH Ti' |43| sample which was used in the Monte 
Carlo simulations |27l | . It is a typical Lunar Basalt and is 
similar to samples 70017, 70035 in ref. For the 'A-15 
Pyroxene' |^ sample we determined = 23.2 g/cm^ 
and Xq = 24.7 g/cm^ for some Apollo-16 samples fSsj . 
In order to obtain an estimate of the sensitivity of the 
acceptance to the stopping power we also show in Fig. [7| 
the results of a calculation where the angular spread is 
reduced by a factor, equal to the reduction of the shower 
length, (24.7/22.1). Since the detection limits arc pro- 
portional to the third power of the angular spread, this 
results in only a 40% change. 

An additional source of model dependence lies in prin- 
ciple in the details of modelling the roughness of the lu- 
nar surface. For the GLUE experiment, at much shorter 
wavelength, this roughness gave rise to a considerable 
broadening of the angular acceptance and thus to a large 
increase in the acceptance. Since, for the wavelength of 
interest for the LOFAR and WRST telescopes, the angu- 
lar spread is already large, the additional effects of sur- 
face roughness can be ignored as they give rise to only a 
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minor increase in the detection probability. 

The Puma-II back-end which is available at the WRST 
telescope allows for storing all data on disk for each ob- 
servation period. A total of 12 hours of observing time 
can be stored. The analysis of the data can thus be done 
off-line allowing full flexibility in optimizing the correc- 
tions for ionospheric dispersion. Because of the Nyquist 
sampling of the signal in the PUMA-II back-end the pulse 
may cover -after correction for dispersion- two to three 
sampling times. In the off-line analysis this can be taken 
into account. In the LOFAR operation the signal is also 
stored but only after a trigger condition is met in a fast 
on-line signal analysis. 



VI. SUMMARY 



hundreds or even thousands of km, making the moon an 
excellent detector. The difference in the time of arrival 
of the two particles determines the mass of the original 
cosmic-ray nucleus. 
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We have demonstrated the clear advantage of using 
radio waves at frequencies well below the Cerenkov max- 
imum. The optimum frequency will be that where the 
length of the shower, of the order of several meters in the 
lunar regolith, is of the same order of magnitude as the 
wavelength of the radio waves where the radio-emission 
pattern is nearly isotropic. The advantage of going to 
lower frequencies applies to all experiments where the 
radiation crosses a boundary between a dense medium 
to one with a considerably lower index of refraction. 

We have shown that the gain in efficiency at lower fre- 
quencies is such that with the upcoming LOPAR facility 
one can seriously investigate realistic top-down scenarios 
for UHP neutrinos and be sensitive to neutrino fluxes well 
below the Waxman-Bahcall limit. Pven now with the 
existing WSRT, profiting from its capability to measure 
right in the radio-frequency window where the detection 
efficiency is highest, one is able to set limits on neutrino 
fluxes orders of magnitude below the present limit in only 
a 100 h observation period. 

Por UHE cosmic rays the LOPAR facility offers, be- 
cause of the availability of an optimal radio-frequency 
window, a very powerful tool to determine the flux be- 
yond the GZK limit. In only a 30 day observing period 
one is sensitive to a flux which is more than one order of 
magnitude below the extrapolation of the measured flux 
from below the GZK limit. Since the beam of LOPAR 
is determined by software, much longer observation pe- 
riods should also be attainable, resulting in a sensitivity 
to even lower fluxes. Assuming the GZK limit is real, 
this offers the exciting possibility to measure the density 
of sources for UHE cosmic rays within a range of about 
10 Mpc. 

As an additional topic one may access the composition 
of cosmic-rays (proton versus heavy nuclei) by search- 
ing for the predicted coulomb dissociation of heavy nu- 
clei passing in the neighborhood of the high density pho- 
ton fleld of the Sun, the so called Gerasimova-Zatsepin 
effect l38l . The original predictions were recently re- 
vised [33 ; showing that the separation of the two dissoci- 
ated daughter nuclei at 1 AU from the Sun is as large as 



APPENDIX A: ANGULAR SPREAD 

Since the angular spread of the intensity of the 
Cerenkov radiation around the Cerenkov angle for the 
case X fa L is crucial for our considerations, we will 
present here a discussion of this case which is indepen- 
dent of the parameterizations given in the literature at 
shorter wavelength. 

In the literature the intensity of Cerenkov radiation 
from a hadronic shower, with energy E^, in the lunar 
regolith, in a bandwidth Ai/ at a frequency v and an 
angle 9, has been parameterized, based on Monte Carlo 
simulations, as [llLlla .llq 

F{e,u,E,) ^ 3.86X 104e-((''-''=)/^^)'(^^)' 

"^11020 eV/ + 
"(10^) 

where vq = 2.5 GHz and the spreading Ac is given by 
Eq. 101. Eq. IjAip has been experimentally shown ^3 to 
be accurate within a factor 2 for a shower of particles that 
would correspond to a primary particle of about 10^^ eV 
at frequencies exceeding 1 GHz where the wavelength is 
small compared to the longitudinal extent (length) of the 
shower. Por the case in which the wavelength is compara- 
ble to the length of the shower, of interest for the present 
investigation, one may expect deviations from the simple 
parametrization where our main concern is the angular 
spread of the Cerenkov radiation. 

To focus on the angular spread, we have derived the 
angular distributions for two different shower profiles fol- 
lowing the approach given in For the first one, called 
the "block" profile, the number of charged particles is 
constant over the shower length L ^ Lb, Pb{x) = 1 for 
Q < X < Lb. This profile is not realistic for a shower 
as the full intensity suddenly appears and disappears. 
Por this reason we have also investigated a second pro- 
file where the charge in the shower appears and disap- 
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FIG. 9: The angular spread around the Cerenkov angle for 
the block and the sine shower-profile functions are compared 
to the parametrization used in this work (e: Eq. IIAIL g: 
Eq. igSj; b: Eq. (glj; s: Eq. JXH). The left (right) hand 
displays the results for 2.2 GHz (lOOMHz) respectively. 



pears following a sine profile, Ps{x) = sinirx/Ls with 
< X < Ls. 

For the block longitudinal profile we reproduce the well 
known intensity distribution found by Tamm |47l | for a 
finite length shower, normalized to unity at the Cerenkov 
angle, 



h{0) = 



sni a sm irx 



TTX 



with 



X = (cos 6* - l/n)L/X 



For the normalized "sine" profile we obtain 



sinfl 



cos TTX 



(l-2x)(l + 2x) 



(A2) 



(A3) 



(A4) 



The predictions of these two formulas are compared with 
the parametrization of Eq. at 2.2 GHz and 100 MHz 
for a shower of 10^° eV. To reproduce the angular spread 
of this calculation at 2.2 GHz we choose Li, = 2.5 m and 
= Lf, X 4/3 = 3.4 m, see Fig. El The results are 
also compared to those of the gaussian parametrization 
proposed in [T^ . 



sinO 



smf 



(A5) 



with Z given by Eq. (|A7|) . The value for Zq, Eq. 
is chosen such that for a small angle expansion around 
0c it agrees exactly with Eq. (|A1|) . From the figure it is 



seen that the simple parametrization of Eq. (|Aip is repro- 
duced well by all three analytic forms. Eq. (|A2() shows 
the well known secondary interference maxima, due to 
the sharp edges of the profile , which are not realistic 
for our case. Keeping parameters fixed the angular distri- 
butions are now compared at 100 MHz (right hand panel 
of Fig. O. The three analytic forms. Equations ljA2p . 
(jA4p . and IjASp agree quite accurately but differ consid- 
erably from Eq. IjAip . The reason for this difference lies 
mainly in the pre- factor sin^ 9, missing in Eq. (|A1() . which 
accounts for the radiation being polarized parallel to the 
shower and thus that emission at 0° and 180° is not pos- 
sible. 

On the basis of the arguments given above we will use 
the gaussian parametrization, accurate at small and large 
angles, 

F{e, u,E,)= 3.86 X 10^ e-^' f ^ V f ' (A6) 

Vsmt^c^ \ a J 



1020 eV/ \vo{l + {v/voY-'^'^)J '100 MHz 



1.44-j) (" 



with 



Z = (cos 6* - l/n)Zo 



The value for 



180 



Jy, 



(A7) 



(A8) 



whth Ac measured in degrees (see Eq. For small 

spreading angles around the Cerenkov angle, i.e. short 
wavelengths, this expression agrees with the results of 
Monte Carlo simulations 0, 0, 0| while at large wave 
lengths the formula agrees with the analytic results [i^ . 

As a last point we compare the length of the shower, 
1.7 m according to Eq. with the value for the length 
used in Eq. (|A4|) . Ls = 3.4 m. It should be reahzed 
that for the sine profile, only for half its length (i.e. 1.7 
m) the density of charged particles exceeds 70% of the 
maximum value, which is the definition of the shower 
length in the Monte-Carlo simulations. The agreement 
is thus excellent. 



APPENDIX B: SHALLOW SHOWERS 

In this paper we have implicitly assumed that the 
showers develop well inside the lunar regolith. For the 
emission of long wavelength radio waves from cosmic ray 
showers which are close to the lunar surface this assump- 
tion needs a more detailed consideration. The proximity 
of the surface, through mechanisms like mirror charges, 
could severely diminish the amount of Cerenkov radiation 
through the surface. 

Even though the general problem has not been studied, 
the -in some sense- inverse problem has been studied |46j , 
namely that of an electron beam in close proximity to a 
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dielectric. When the velocity of the electrons exceeds 
that of the velocity of light in the medium (cm = c/n, 
where n is the index of refraction) Cerenkov radiation 
is induced in the medium. The occurrence of this pro- 
cess has been verified experimentally ^3 ■ In Rsf. the 
amount of Cerenkov radiation is calculated as a function 
of all key parameters in the problem such as the distance 
a of the electron beam from the surface, the electron ve- 
locity /3 = v/c, the wavelength of the radiation A, and 
the angle 77 of the radiation with respect to the surface of 
the dielectric. In the limit of ultra-relativistic particles 
(/3 = 1) the dependence of the intensity on a reads 



W 



-4ir aVn^ — 1/A 



(Bl) 



This equation differs from that quoted in |4l| where in- 
stead the limit /? ^ 1/n has been used. Eq. ljBl(l clearly 
shows that proximity of the surface is only an issue when 
AnaVn^ - 1/A < 1 or for a < X/iAny/n^ ~ 1) = 0.05A 
for the lunar rcgolith. This implies that this effect should 
be considered only for showers making very small an- 
gles to the surface, i.e. "shallow showers" with sin9 < 
0.05A/(0.5 L), where L is the shower length. For L ~ X 
this corresponds to an angle of less than 6° . 

To test the contribution of shallow showers to our re- 
sults wc have made calculations in which the contribution 
of shallow showers are excluded. The difference in the re- 
sults is barely visible in the plots and this effect can thus 
safely be ignored. 

It should however be realized that it is only an extrap- 
olation to apply the conclusions based on the work of 
Ref.[46l| to the present problem. A thorough theoretical 
treatment should be performed. 



APPENDIX C: LUNAR REGOLITH 

The properties of the regolith play an important role 
in the present calculations. The index of refraction and 
the loss tangent (determining the attenuation length of 
radio waves) have been determined from samples brought 
back to Earth in the Apollo missions. From the average 
values one extracts an index of refraction of n = 1.8 and 
an attenuation length of A^ = (9/i/[GIIz]) m [MEl for 
the intensity of radio waves in regolith. 

One issue of particular interest for the detection ef- 
ficiency of neutrino-induced showers is the thickness of 
the regolith. In his thesis Takahashi made an ex- 
tensive study of the absorbtion and reflection of radio 
waves in the frequency domain of 0.1 to 10 MHz for real- 
istic depth profiles based on • expects a smoothly 
varying attenuation length of A,- = 4.08(i^[GHz])~ ®^ m 
for depths ranging from 200 m up to 100 km. At these 
depths one expects a decreased attenuation length due 
to the increased density of the material. Since the den- 
sity of rock is about twice as high as that of the regolith 
also the value of the index of refraction (and thus the 
Cerenkov angle) and the length of hadronic shower (and 
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FIG. 10: The dependence of the detection limits on properties 
of the lunar regolith. The solid curve represents a calculation 
made as if the regolith extends to a depth of 500 m, with only 
20 m assumed for the dotted curve, while pure lunar rock is 
taken for the short dashed calculation. 



thus the spread around the Cerenkov angle) should be 
modified. To study these different effects we compare in 
Fig- El the results of three diS'erent calculations. In the 
first (the solid curve) we assume that the properties of 
the rcgolith layer apply to depths of 500 m. In the sec- 
ond (dotted curve) the rcgolith expends to a depth of 20 
m and we assume that no radiation from deeper layers 
reaches the surface. In the third calculation we assume 
that the lunar rock extends to the surface of the moon. 
The density of the rock is taken twice that of the regolith, 
the index of refraction equal to 2.6, the radio absorbtion 
length equal to one third of that of the rcgolith, which 
takes into account the larger value of the loss tangent in 
rock. 

At a frequency of 2.2 GHz one finds, by comparing 
the drawn and the dashed curves, that only showers in 
the upper part of the regolith are detectable due to the 
relatively strong absorbtion of radio waves. At this fre- 
quency the detectability of showers in rock is much higher 
due to the fact that the spread around the Cerenkov an- 
gle is twice as large due to the reduced length of the 
shower which in turn is a consequence of the larger den- 
sity. However the rock is at most areas covered by a 
layer of regolith which will completely attenuate the ra- 
dio waves. The calculations for 20 or 500 m of regolith 
(hardly any difference between the two) thus yield a con- 
servative lower limit. At 100 MHz the situation is quite 
different. The calculations for pure rock or pure regolith 
give rise to very similar limits. The increase in spreading 
width for the calculation in rock is apparently compen- 
sated by an increased attenuation. Due to the larger 
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wavelength the contribution of showers deep under the 
surface are important. If only the showers in the upper 
20 m of regolith are taken into account in the calculation, 
the limit changes by about 1 order of magnitude. This 
calculation completely ignores the emission from deeper 
showers which is clearly unrealistic. If one adds the con- 
tribution from the deeper rock layer one will obtain a re- 
sult close to the drawn curve. Including the fact that for 
rock covered by a layer of regolith the index of refraction 
changes more gradually, giving rise to reduced reflection 



near the surface, would even give a more stringent lower 
limit. Again, the calculation taking a maximum depth of 
500 m gives a realistic estimate for the limit. 

The calculations in the main part of this work account 
for radiation emitted from a depth of not more that 500 
m. On the basis of the arguments presented in the above 
it will be clear that this gives a realistic estimate for lower 
bounds on neutrino fluxes. For cosmic rays none of these 
considerations are important as all induced showers lie in 
the upper part of the regolith. 
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